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SUMMARY 


A Doppler tracking method is presented in which a very small , 
simple CW beacon transmitter is used with four Doppler receiving sta- 
tions to obtain the position and velocity of a space research vehicle. 
The exact transmitter frequency need not be known, but an initial posi- 
tion is required, and Doppler frequencies must be measured with extreme 
accuracy. The errors in the system are analyzed and general formulas 
are derived for position and velocity errors. The proper location of 
receiving stations is discussed, a rule for avoiding infinite errors is 
given, and error charts for ideal station configurations are presented. 
The effect of the index of refraction is also investigated. The system 
is capable of determining transmitter position within 1,000 feet at a 
range of 200 miles. 


INTRODUCTION 


The accurate determination of the position and velocity of a small 
research vehicle is an important problem. Many elaborate methods for 
making this determination have been developed, most of which require 
either high-powered, complicated, and expensive ground equipment or 
complicated airborne equipment (ref. l). There is need for a tracking 
system which uses extremely small simple equipment in the vehicle and 
moderately simple inexpensive equipment on the ground. Such a system 
is afforded by transmitting a low- power CW signal from a radio beacon 
in the vehicle and receiving Doppler signals at several ground stations. 

The simplest Doppler system of this type requires a very stable 
transmitter and three receiving stations. The disadvantage of this 
method is the precision with which the transmitter frequency must be 
known. By using four stations, the problem becomes more involved, but 
it is possible to find the trajectory even though the transmitter fre- 
quency is not known precisely or even if it shifts. In this investi- 
gation the problem of finding the position by using four stations is 
solved. Position- error formulas are derived and optimum station 
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locations to minimize these errors are discussed. The accuracies 
required in measuring Doppler frequencies and the effect of the index 
of refraction are discussed. Also, the velocity and errors in velocity 
are found. 


SYMBOLS 


A, B, C Doppler receiving stations 

A x;B x ,C x ;Ay, 


By, Cy - . . 


a, b, c 
a,b, c 
a x c xJ a y* 

by,CyJ . . . 

aa 


cofactors of the elements a x } b x , c x ; ay, by, 

Cy; . . * respectively, in determinant E 

distance from stations A, B, and C, respectively, to origin 
vector position of stations A, B, and C, respectively 


D = 


E = 


5fa 

d g b 

Ss c. 

8s a 

dx 9 

5x 9 

8x ' 

dy ' 


Ss-l. 

- — , - — ; . . . (see eq. ( 28 )) 

oy oy 


probable error in x a , Xb, or x c ; y a , yb> or Vc} 
when these probable errors are all equal 

velocity of propagation of radio waves in free space 
velocity of propagation at source 


x a 

y a 

s a 

x b 

^b 

s b 

x c 

y c 

s c 



1 

a z 

b x 


b z 

c x 

c y 

c z 


(see eq. (A2)) 


vapor pressure, mb 
transmitter frequency 
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f a ,f b ,f c ,f 0 received, frequencies at stations A, B, C, and 0, 
respectively 

Af error in f 

of probable error in (f Q - f a ), ( f o “ ^b); or (^o " f c) ^ en 

these probable errors are all equal 

d ( f o - f a)> d ( f o - f b)> d ( f o - f c) differential of (or error in) 

(fo _ ^a)> ( f o " ^b)-> and (*o “ ^c) 

quantities defined in eqs. (l6), (l7)> and (l8), respectively 

length of sides in rectangular configuration 
electron concentration, per cm^ 
index of refraction at transmitter 
station 0 at the origin 
atmospheric pressure, mb 

distance from origin to transmitter or distance from sta- 
tion 0 to transmitter 

vector position of transmitter relative to origin 
differential of vector position of transmitter 
probable error in position of transmitter 

distance from stations A, B, and C, respectively, to trans- 
mitter (see eqs. (l0), (ll), and (12)) 

vector position of transmitter relative to stations A, B, 
and C, respectively 

r ao , ^o* r co* r oo initial distances from stations A, B, C, and 0, 

respectively, to transmitter at time t Q 

Ar a error in r a 

S location of source (transmitter) 


Z,m 

N 

n 

0 

P 

r S Tr, 


r 


dr 


ar 


r a' r b' r c 


r a > ^ r c 


k 


S c 


S ■ Si i s 

a 9 b* c 


ds a' ds V ds c 


as 

as a , os^y os Q 

T 

t 


v a> v b> v c 


v x^V V Z 


ov 


crv , crv , crv 
x 7 y z 


X a> X a^ X b' X bJ 

Wc 


dx, dy, dz 
crx, cry, o z 


cofactors of s a , s-^, and s c , respectively, in 
determinant D 

range differences, r a - r (eq. (23)), r b - r, and r Q - r, 
respectively 

differentials of (or errors in) s a , s^, and s c , 
respectively 

probable error in s ft , s^, or s c when these probable 
errors are all equal 

probable error in s a , s^, and s Q , respectively 

temperature, °K 
time 

initial time 

radial components of velocity of transmitter (source) rela- 
tive to stations A, B, and C, respectively 

x, y, and z components, respectively, of transmitter velocity 

probable error in transmitter velocity 

probable errors in v x , Vy, and v z , respectively 

cofactors of x a , y a ; x-^, y^; x c , y c in determinant D 

coordinates of transmitter 

differentials of (or error in) x, y, and z, respectively 
probable error in x, y, and z, respectively 


x b'yb'*b * ■ * 


coordinates of stations A, B, . . . 


differentials of (or errors in) 
x b> yb> z b^ * * * 


"a' 




x a> 


-T'K-h ro h t- 1 
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ox a> oy a , Gz a ; probable errors in x a , y a , z a ; y b , z b ; . . 

ox b ,ay v az b> • • • 

5 propagation error factor (see eq. (4j)) 

e error factor, ratio of probable error in r to probable 

error in range differences, ar/crs 

t| per unit error in index of refraction 

A wavelength of transmitter signal at source 


PRELIMINARY ANALYSIS 


A brief analysis of the three- station method is presented. This 
analysis shows that a range measurement depends on the transmitter fre- 
quency. This frequency dependence can be eliminated by measuring range 
differences between stations. Since at least three range differences 
or four stations are required, an analysis of the four- station method 
is presented. 


Three- Station Method 

A nonrelativistic approximation of the Doppler equation for the 
frequency observed by a receiver at a fixed point A is given by 



where v Q is the radial 

d 

the point A, and c s is 
Thus, 


or 


velocity component of the source relative to 
the velocity of propagation at the source. 


- T ( f - 

(2) 

= f a) 

(3) 


v a = K f - f a) 


(4) 
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If the initial position of the source is known, the subsequent 
radial distance (range) from point A can be determined from 


r 


a 


ao 



A(f - f a )dt 


( 5 ) 


With two receiving stations A and B, the ranges r a and r^ are 

then known. Figure 1 illustrates that the intersection of two spheres 
with radii r a and r^ locate the transmitter somewhere along a circle; 
hence, an additional station C (located off the axis A, B) is required 
to locate the transmitter. (The intersection of the sphere C with the 
circle actually gives two positions, but in a practical case one of 
these may be disregarded since it is beneath the surface of the earth. ) 

In practice the received frequencies can be measured very accurately, 
but the transmitter frequency f cannot be known exactly because it may 
shift during flight. It is difficult to keep the frequency of a very 
small transmitter constant under the extreme environment encountered in 
a rocket- propelled vehicle. From equation ( 5 ) it can be seen that an 
error Af in f gives an error in r a of 


L 

1 




A Af dt 


( 6 ) 


As an illustration of the magnitudes involved, consider a typical 
case where f = 200 megacycles (A ~ 5 feet) with a constant error of 
100 cycles per second. Over a time interval of 300 seconds the error 
in r a would be 150,000 feet. 


Four- Station Method 

The error due to transmitter- frequency change may be eliminated by 
measuring the frequency difference between two Doppler receiving sta- 
tions. Let 0 represent an additional receiving station located at 
the origin; the range difference is 


= r 


ao 


- r. 


00 



(T) 


C\J 


ro H 


7 


There will be a small error in A since f is not known; but f 
differs from f Q only by the amount of the Doppler shift which is less 

than 5 kilocycles for a transmitter speed of 25,000 feet/sec. Thus, 
if A is based on f G the error in s a is less than 0.0025 percent 

(or less than 26 feet at a range difference of 200 miles). 

There may also be a small error in s a because of the difference 

in time it takes for the signal to reach different receivers. If the 
receiving stations are 200 miles apart, the maximum time difference is 
only 1 millisecond. If the velocity of the beacon suddenly changed by 
25,000 ft/sec, which gives a 5-kilocycle change in Doppler signal, the 
maximum error in s a would be only 25 feet. (Both these errors can 

be eliminated, if desired, once the position and velocity of the trans- 
mitter are determined. ) Thus, by measuring the difference in frequency 
between two stations, the range differences can be determined without 
knowing the transmitter frequency. 

If it is assumed that the range differences s a , s b , and s c can 
be measured accurately, it can be seen that four receiving stations are 
required to locate the transmitter. In figure 2 receiving stations are 
shown at A, B, C, and 0, with the source at S. A sphere with radius r Q 

is illustrated with S as the center. Since the distance from S 
to A is r a , then 

r a = ( r a - r o) + r o 


Hence, the sphere with radius s a and center A is shown to be tangent 
to the sphere with center S. Similarly, spheres with radii s b and 
s c are shown around centers B and C. Now, given the range differ- 
ences s a , s b , and s c , S is located at the center of the sphere 
which is tangent to the spheres A, B, and C and goes through the point 0 
A geometric solution is thus obtained. 


FOUR- STATION SOLUTION 


The geometrical solution just presented is of little value in 
obtaining the actual solution. This solution may be obtained by writing 
the equations for the four ranges (fig. 5) as 


(9) 
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r a 2 = ( x - x a) 2 + (y - y a ) 2 + ( z - z a ) 2 (10) 

r b 2 = ( x - x b) 2 + (y - y b ) 2 + ( z - z b) 2 ( 11 ) 

r c 2 = ( x - x c) 2 + (y - y c ) 2 + ( z - z c) 2 (i2) 

Subtracting equation (9) from equation (10) yields 

r a 2 - r 2 = -2x a x - 2y a y - 2z a z + x & 2 + y& 2 + z& 2 (l 3 ) 

With the substitutions 

2 2 p p 

a = x a + y a + z a (14) 


r a 2 - r2 = (r a + r)(r a - r) 

= ( 2r + s a) s a 
= s a 2 + 2rs a 


equation (lj) 


becomes 

x a x + y a y + s a r 


Similarly, 


x b x + yby + s b r 

x c x + y c y + s c r 


a 2 - s a 2 - 2z a z 
2 

2 p - 

b - s b - 2z b z 

2 

? P 

c - s c - 2z c z 
2 


= *a 
= K b 


( 15 ) 


(16) 


( 17 ) 

(18) 


In equations (l6) to (l8) a, x a , y a , z a ; b, x b , . . . are 
constants which depend on station locations; s a , s-j^, and s c are 
known data from the Doppler stations; x, y, and r are unknowns to 
be solved. Actually z is also an unknown, but since z a , z b , and z c 

are usually small, an approximate value may be assumed in order to obtain 
a first approximation of x, y, and r. 


ro H f 
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The solution is then (using cofactors of determinant D) 


K a X a + K b X b + K C X C 

x = 

D 

(19) 

K a Y a + K b Y b + KcY c 
y D 

(20) 

KaSa + K b S b + KcS c 

r = 

D 

(21) 

z = - y^ 

(22) 

The computed value of z may be compared with the assumed value, and 
if it is different it may be inserted in equations (l6), (17), and (l8) 
to find more exact values of x, y, and z. The process converges 
rapidly. In several actual cases where the first assumed value of z 
was far from the actual value, the process converged in fewer than six 
iterations. In a typical case in which z was actually 500,000 feet 
but was assumed to be zero, the errors in z for successive iterations 
were 32, 200, 2,100, 132, and 8 feet. 


ERRORS IN POSITION 


This four- station solution does not always give the transmitter 
position. For example, consider four stations located at the corners 
of a square, with the transmitter along a perpendicular through the 
center of the square. No matter where the transmitter is located along 
this perpendicular, r = r a = = r c and, hence, s a = s-k = s c = 0; 

thus, the position is indeterminate and large errors occur in the vicin- 
ity of the perpendicular. The station configuration has a pronounced 
effect on the errors. 

In order to investigate the errors it will he convenient to derive 
an error factor e defined as the ratio of the probable error in the 
position determination to the probable error in the range- difference 
measurements. The error factor will be examined for conditions which 
cause infinite error, and error- factor contour plots of transmitter 
position will be given for various station configurations. The posi- 
tion error depends on range- difference measurements (which in turn 
depend on frequency measurements and propagation effects ) and on measure- 
ments of station positions. All these errors will be discussed herein. 
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Errors Due to Range- Difference Measurements 

Derivation of error factor .- Of primary importance is the error in 
position due to errors in the measurement of the three range differences 
s , s^, and s c . The quantity 


s 


a 


r 


= f - x a ) 2 + (y - y a ) 2 + (z - z a ) 2 - ^X 2 + y 2 + z 2 (23) 

depends on the transmitter coordinates and the station locations. 

Consider the variation ds a in s a caused by the variations dx, 
dy, and dz. Thus, 


dso 8s n 8s a 

ds n = — - dx + — - dy + — - dz 
a 8x Sy 8z 

For simplicity this equation may be written in the form 

ds a = a x dx + ay dy + a z dz 


(24) 


(25) 


Similarly, 


and 


ds b = b x dx + by dy + b z dz 


ds c = c x dx + Cy dy + c z dz 


(26) 

(27) 


where 


a x 


d g a 

8x 

x - *a x 
r a " r 



(28) 


With the other coefficients similar in form. 


r— | CVJ 



ro i — 1 
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Then, the deviations in x, y, and z due to ds a , ds b , and ds f 
are given by solving equations ( 25 ), (26), and (2 7 ) as follows: 


‘ T ds a + T ds f + F ds = 


, A y , By C v 

dy = E dSa + ^ ds b + ^ ds 


E 


dZ = E dS a + F dSd + T dSc 


(29) 

(50) 

(5D 


where A x , By, and so forth are the cofactors of a x , b y , and so forth 
m the determinant E. J 


The probable errors in the range differences may be expected to h P 
u£ S* uncOTr 5 ate ‘ i because each station is similar but^ndipendenj 


crs = as = as-u = as. 


(52) 


The probable errors in the position coordinates 


are then given by 


II 

OJ 

1 

(Ax 2 + B * 2 + 

OJ 

o* 

) (crs) 2 
E 2 

(55) 

(cry) 2 = 

('V 2 + ®y 2 + 

c y 2 ) 

i (os) 2 
E 2 

(3*0 

(az) 2 = 

(Az 2 + B z 2 + 

Ct 2 ) 

(as) 2 

E 2 

(55) 


NoWj since 


— > 2 p p o 

dr = dx + dy + dz^ 


the probable error in the vector r is given by 

ar = \J(ox) 2 + (ay) 2 + (az) 2 


(56) 


(57) 
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Thus, 


€ 


or 

os 


= Ha- 

E' ' 


x 2 + 


B x 2 + C x 2 + 


p p 

A/ + By 


+ Cy + A Z + 


B z 2 + 


( 58 ) 


This error factor depends only on the station configuration and 
the transmitter position. 

L 

Infinite errors - triangle rule .- Infinite errors occur whenever E, ^ 

the demoninator of the error factor, is zero. For the special case when 2 

all the receivers are in the x-y plane (this is the only case which will j 

be considered) appendix A shows that infinite error occurs whenever the 4 

transmitter is in such a position that either z = 0 (i.e., in the 

x-y plane) or D = 0 . 

In appendix B and figure k the station configurations which allow 
D = 0 are studied; from this study the following triangle rule of infi- 
nite error is established: Infinite error cannot occur (except in the 

x-y plane) if any one station lies within or on the triangle formed by 
lines joining the other three stations (fig. 5 ); if no station lies 
within such a triangle, then infinite error exists. 

As an illustration, the special case of a rectangular configuration 
(fig. 6 ) is considered in appendix C, and it is shown that infinite error 
occurs on the two vertical planes which are perpendicular bisectors of 
the sides of the rectangle. For other more intricate configurations, 
the determination of the locations in which infinite error occurs becomes 
much more complex. 

Error factors for various configurations .- In addition to knowing 
the existence of infinite errors, it is important to know the actual 
magnitude of the errors. It may be surmised from the infinite- error 
analysis that a desirable station configuration would occur when one 
station is at the center of an equilateral triangle formed by the other 
three. This equilateral- triangle configuration was analyzed by computing 
the error factor from equation ( 38 ). (The receiving stations have been 
assumed to be located on the earth 1 s surface rather than on a flat plane.) 

The results are shown in figure 7 in the form of contour plots of con- 
stant error factor at various altitudes. This configuration is an ideal 
one for general applications where the azimuth angle may not be known 
in advance. 

Another important configuration is the isosceles triangle with one 
station at the midpoint of the base. Figure 8 shows contour plots for 
this case (the earth 1 s curvature is again taken into account). It may be 
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noticed that, in general, the overall error factor has not been improved. 
However, in one direction (over the apex of the triangle) the error fac- 
tor has been improved at the expense of making the error factor in other 
directions greater. This configuration is useful when the approximate 
azimuth direction is known in advance. A check on the error- factor 
formula (eq. ( 38 )) was made by substituting (in eqs. ( 19 ) to (22)) values 
for s a , s-^, and s c which were in error by 100 feet. This check was 

made for the isosceles- triangle configuration at an altitude of 
500,000 feet. The results checked within 0.2 percent in the range 
covered along the axis 20° counterclockwise from station A (fig. 8(c)). 

One particular configuration for use at NASA Wallops Station 
required land stations which were readily accessible. The stations were 
placed at Wallops Island, Va.; Cape Hatteras, N.C.; Langley Field, Va.; 
and Dover, Del. Error contour plots for this case are shown in figure 
One azimuth direction, as expected, gives infinite error, but in the 
direction of the predicted transmitter position (1^5° clockwise from 
north), the error factor is acceptable. 

Frequency accuracy .- An examination of the contour plots indicates 
that an error factor of about 10 can be achieved. If position data 
within 1,000 feet are required, as must be within 100 feet. This 
accuracy requires the probable error in the total number of cycles 
counted to be within 20 cycles (X = 5)- If this count were over a 
period of 300 seconds, the frequency accuracy required would 0.07 cycle 

per second, an accuracy of 1 part in 3 X 10^. 

It is feasible to use local oscillators which have stabilities of 

1 part in 10^ drift per day at each station. If the oscillators are 
checked against each other within 8 hours, the proper accuracy could 
be obtained. A desirable method of calibration is to check all stations 
simultaneously against a common source (whose frequency need not be 
known exactly). 

Propagation errors .- The index of refraction may vary and cause 
errors. In the troposphere the index changes because of air density 
variations according to (ref. 2) 


77.6 / 48l0e\ c , . 

n = 1 + -T“ ( p + ~T~) 10 (59) 

From the ground up to 50,000 feet, the index varies from about 1.0005 
to 1.0000. The index may be calculated from radiosonde data of pressure, 
temperature, and humidity as a function of altitude. 
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In the ionosphere the index varies because of ionization and 
depends on the electron density and frequency (ref. 3) as follows 
(where f is measured in kilocycles): 


r sin / * . 

n . - — (W) 

The electron density becomes important at an altitude of about 
200,000 feet; at 1 x 10^ feet it attains a maximum density of the order 
of 1.5 x 10^ per cubic centimeter (ref. 4). Thus, at 200 megacycles 

the index varies from 1.0000 to 0.9985 between 200,000 and 1 x 10^ feet. 

It appears that a frequency of at least 400 megacycles should be used, 
in which case the variation would be only 0.0004. The index may be cal- 
culated by using electron density measurements from rocket soundings. 

The effect of the Index of refraction on position determination 
can be seen from equation ( 7 ) where A = c Q jnf. However, the error does 

not enter in the same manner as the frequency error in ^f 0 - f a ) because 
the errors in s a , s-^, and s c due to errors in n are correlated. 

If the simplifying assumption is made that the per unit error in n is 
a constant q, and if the initial range difference is neglected, then 
the per unit error in s a , s^, and s c is rj and the error in x is 

(from eq. ( 29 )) as follows: 

Ajf Ryf Cy 

dx = — T)s a + — n s b + -g- O 1 ) 

The probable error in position is then 

tfr = E I (^ xS a "®" ^x s b "** ^x s c) t (4ly-s a + ^y s b t ^y s c) + (-^z s a“^ ^z s b* ^z s c) (42) 


The propagation- error factor may be defined as follows: 


5 = 


ar 


TJ X10 6 


t q— b I " 2 f 

— — y( A x s a + ®x s b +C x s c) + ( A y s a + By s b + ^y s c) + ( A z s a+ B z s b + Cz s c) 


(*o) 


H OJ 



hi H OJ rr\-d- 
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Contour plots of this error factor are shown in figure 10 for the 
equilateral- triangle configuration. If the error in n were a constant 
value of 0.0001, and an error factor 5 of 7 could be achieved, then 
ar would equal 700 feet. 

The variation of the index of refraction with altitude can also 
cause propagation in a curved path, which will produce errors. This 
effect is especially important at low elevation angles; however, it is 
not investigated herein. 


Station- Location Errors 

Another aspect of the problem is the accuracy of measuring the 
location of the receiving stations. Appendix D shows that the ratio of 
probable error in missile position to the probable error in station 
location is equal to the error factor e or 

qr 
qa 

This result may also be seen quite simply from an examination of fig- 
ure 3- Suppose station A were moved a small distance. It is seen that 
a change in the vector a* (vector position of A) is the same as a 
change in the vector r a . Hence, a change in a* must have the same 

effect as a change in r a and, consequently, in s a - 

If the error in station location were 5 feet and if the error fac- 
tor were as great as 5 0 , the error in transmitter location would be 
only 250 feet. 



VELOCITY DETERMINATIONS 


The velocity of the transmitter may be obtained in two different 
ways. In the most direct method the velocity components are obtained 
from the x, y, and z positions as a function time. In another method 
the Doppler frequencies are used in conjunction with the transmitter 
location. Appendix E shows that the velocity components are 

v x = ^ ~ + ®x( f ° ” ^b) + ^x(^o “ fc)] (^5) 

v y - | [%( f 0 - fa) + By(f 0 - f b ) + Cy(f 0 - f 0 )] 


( 46 ) 
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A 

E 



fa) + (f o - f*b) + ^z(fo ■* fc) 


( 47 ) 


where A x , By, and so forth depend on receiving- station configuration 

and transmitter location. Actually, the two methods differ only by the 
errors introduced in the data- reduction process and, hence, preference 
for one method will be determined by this process. 


ERRORS IN VELOCITY L 

1 

2 

If it is assumed that the position error is negligible, then the 3 

probable error in velocity (see appendix F) is given by 4 

ov = Aeaf (48) 

The instantaneous frequency could not be expected to be as accurate 
as the integrated frequency. If an accuracy of 0.5 cycle per second 
can be obtained, then velocity data would be within (5)(l0)(0.5) = 25 feet 
per second. 


CONCLUSIONS 


An analysis has shown that a four-station Doppler system is capable 
of giving position data on a research vehicle to within 1,000 feet 
(except at low elevation angles) by using a simple radio beacon trans- 
mitter whose frequency need not be known accurately. However, Doppler 
frequencies must be measured very accurately and variations of the index 
of refraction should be taken into account for best accuracy. 

Errors for the equilateral- and isosceles- triangle configurations 
may be found from the contour plots given; for other configurations, 
formulas may be used. Infinite errors may be avoided by locating one 
station inside the triangle formed by the other three stations . 

The equi lateral- triangle configuration with one station in the center 
is the best configuration for general use, although the isosceles- 
triangle configuration with a station on the midpoint of the base will 
give better results in a specific direction. Other configurations may 
be satisfactory but should be checked with the error formulas. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., January 10, 1961. 
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IN. 


APPENDIX A 

EXISTENCE OF INFINITE ERRORS 


If all stations are in the x-y plane, it will be shown that infinite 
error occurs when D = 0. Infinite error occurs whenever E, in equa- 
tion (38), is zero. The determinant E is given by 


[1 

L 

2 

5 

+ 


a x 

a y 

&Z 


8s a 

bx 

6s a 

8y 

ds a 

bz 

^x 

b y 

bz 

= 

8s b 

bx 

8s b 

8y 

ds b 

dz 

C X 

c y 

c z 


ds c 

bx 

8s c 

8y 

ds c 
bz f 


(Al) 


By using equation (28) 


the following determinant may be 


obtained: 



The next step is to factor x } y, and z out of the columns and r a , 
r^^ and r c out of the rows. For the special case when z a = zb = z c = 0 
(i.e., all stations located in the x-y plane) E becomes 



1 

_ ^ 

X 

_ fa 
r 

1 

_ Za 
y 

_ fa 
r 

1 

. fa 
r 

xyz 

1 

_ fb 

_ fb 

1 

_ Zb 

_ fb 

1 

_ fb 

r a r b r c 

X 

r 


y 

r 


r 


1 

1 

_ 

X 

_ £c 
r 

1 

_ Yc 
y 

. f£ 

r 

1 

_ rc 
r 


(A3) 
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Thus 


Subtraction of the last column from the first two gives 



_ fa 

X 

. Za 
y 

1 

_ £a 
r 

E - xyz - 

. ^ 

. it 

1 

_ fb 

ra^c 

X 

y 

-L 

r 


_ fc 

X 

_ Zc 
y 

1 

_ £c 
r 


(A4) 



x a ya 


- r 


E - z 

x b yt, 

r*b 

— r 

(A5) 

rr a r b r c 



x c yc 

r c 

- r 


E = - 

zD 

rr a r b r c 



(A6) 

error occurs when 

z = 0 or 

D = 

0. 



r\) h tr 1 
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APPENDIX B 

TRIANGLE RULE OF INFINITE ERROR 


In order to prove the triangle rule of infinite error, a general 
triangular configuration of stations A, B, and C lying in a plane (fig. 5) 
is assumed, with the position of station 0, the origin, variable with 
respect to the triangle. From an investigation of the signs of S a , S*,, 

and S c , it will be shown that if 0 lies in the unshaded regions of 
figure 5, the determinant D cannot be zero and infinite error cannot 
exist. When 0 lies in the shaded regions it will be shown that there 
exists a surface of transmitter positions which will cause the determi- 
nant D to equal zero and give infinite error. A generalization of 
figure 5 gives the triangle rule of infinite error which may be stated 
as follows: Infinite error cannot occur if any one station lies within 

(or on) the triangle formed by lines joining the other three stations; 
if no station lies within such a triangle, then infinite errors exist. 

Properties of S a , S-^, and S c 

In order to determine the signs of S a , and S c , consider the 

identity 


x a 

y a 

a 



x b 

yb 

b 

3 0 

(Bl) 



— > 



x c 

y c 

c 




which is readily verified (since stations lie in the x-y plane and 
z a = z b = z c = 0) by subtracting from each element in the last column 

i t im es the corresponding element in the first column and j times 
the corresponding element in the second column, where i and j are 
the unit vectors in the x and y directions, respectively. Thus, 

S a a + S b b* + S c c* = 0 (B2) 


where 
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s a = x by c - x cyt ^ B5) 

Sb = x c^b - x by C W 

s c = x a yb - x b^a ( B 5) 

If S Q = 0, then — = — and the origin 0 lies on BC, or BC 

a ' x b x c 

extended (fig. 5). When S a > 0, point 0 lies on the same side of BC 

as point A, whereas if S a < 0, point 0 lies on the opposite side of BC. 

The signs of St, and S c may be determined in the same manner to 

obtain the remaining information in figure 5 . Note that inside the 

triangle S a , St,* and S c are all positive. On any one side of the 

triangle one of these quantities is zero. 

The identity (B2) may be written in the alternate form 

- — h t = t + S ° ■ ■- (<? - S) (b6) 

Sa + S c S a + S c 


Figure 4 illustrates this identity and shows that if B lies on AC, 


_ Sb = x 

Sa + S c 


or 


Sa + Sq, + S c = 0 


(B7) 


(B8) 


However, if B lies below AC as shown in figures 4 and 5 and if 0 
lies in the unshaded region below B, then 

— > l (B9) 

Sa + Sc 

or, since S a and S c are both negative, 


Sa + Sb + S c > 0 


(BIO) 


ro h t-> 



4=-v*>t ro H f 
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On the other hand, if 0 lies in the shaded region above AC (fig. 5 ) > 

then 


- - fr — <1 

®a + 


and since S a and S c are both positive. 


S a + ^ + S c > 0 


(Bll) 


(B12) 


In a similar manner, it may be shown that equation (B12) is true in all 
the other regions (except on the extended sides of the triangle where 
equation (B8) holds). 


Regions of No Infinite Error 

It was shown in appendix A that for infinite error the following 
condition must exist: 



x a 

y a 

r a 

- r 



D = 

x b 

yb 

r b 

- r 

= 0 

(B15) 


x c 

y c 

r c 

- r 



B = S a (r a - 

r) + 

^(^b 

- r) 

+ S c 

o 

1 

*1 

II 

o 

(Bl4) 


or 


s a r a + Vb + S c r c = ( S a + s b + s c) r 


(B15) 


Identity (B2) may be written 


Sa(? ' ?.) - Sb(? - ? b ) + Se(? - ?=) = 0 


(B16) 
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or 


S a^ a + + S c^c = ( > S a + ^ + S c )r (B17) 

Since S a , S b , and S c are all positive (as they are when station 0 
lies inside the triangle) it follows that 

S a|^a| + Bbftl + S c |r c | ^ (s a + ^ + S c )|?| (Bl8) 

Hence, equation (B15) cannot be satisfied (except when r^, r^, and r c 

are parallel at r = °o), and infinite error cannot exist when station 0 
lies inside the triangle ABC. (This is also true if station 0 lies on 
a side of the triangle.) 

In the unshaded region below station B (fig. 5) equation (B15) may 
be written (in order to emphasize that S a < 0 and S D < 0) 

Vb * ( s a + Sb + s c) r + (' s a) r a + (-Sb) r b (B19) 

and equation (B17) becomes 

^b - ( s a + Sb + S c )r + (-S a )? a + (-Sb)^ (B20) 

Hence , 

Sbl^bl - ( S a + Sb + S c )|r |+ (-S a )|r a | + (-S b )|^| (B2l) 

and no infinite error can exist. Note that in this region station B is 
enclosed by the triangle ACO. In a similar manner it can be shown tha t 
infinite error cannot exist in the other unshaded regions. 


Regions of Infinite Error 

In order to show that infinite error exists when station 0 lies in 
the shaded regions of figure 5> consider D as a function of x, y, and z 
(transmitter position). If D(x,y, z) > 0 at some point in space and 
less than zero at some other point in space, then D must be zero on some 
surface separating these two points (for no matter what path is taken 
between the two points, D must go through zero somewhere along this 
path) and infinite error must exist. The two points which show this 
are the origin 0 and one of the corners of the triangle. Consider the 
shaded region to the left of side AB. Then, at the origin r is zero 
and (since S c < 0) equation (Bl4) becomes 


i—l i — I CM 
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D(0) = S a r a + 

- ( _s c) r c 

(B22) 

and equation (B17) becomes 

^a r a + = ( 

1 

CO 

0 

4 1 
0 ^ 

(B25) 

Thus, 

S a|^a| + = 

(- s =)| ? cl 

(B24) 

and, therefore, 

D(0) > 0 


(B25) 

At the point A, r a is zero and equation (Bl4) becomes 
D(A) = -(S a + ^ + S c )r - (-S c )r c + 

(B26) 

Equation (B17) becomes 

^b r b = (^a + ^b + S c) 

+ (- s cH 

(B27) 

Since all coefficients are positive 

Sb|r b | = ( s a + Sb + s c) 

1*1 + (- S c)| ? cl 

(B28) 


and, therefore, D(A) < 0. Thus, D = 0 somewhere between stations 0 
and A and infinite error exists. The other shaded areas may be analyzed 
in a similar manner. 
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APPENDIX C 

INFINITE ERRORS IN A RECTANGULAR CONFIGURATION 

For a rectangular configuration, consider figure 6. Infinite 
error occurs when 

x a y a r a - r | 0 m r a - r 

D “ x b r b - r = 1 m r b - r 

x c yc r c - r Z 0 r c - r 



r + r b = r a + r c 

(C2) 

or 

r 2 + 2rr b + r b 2 = r a 2 + 2r a r c + r c 2 

(C3) 

But 

2 ? 2 

r a = r c - 2 my + m 

(C4) 


r b 2 = r^ _ 2 lx - 2my + Z 2 + m 2 

(C5) 


r c 2 = r 2 - 2Z x + Z 2 

(c6) 

When 

the 

equations (CJ+), (C5), and (c6) are substituted into 
following equations are obtained: 

equation (C3) 


rr b = r a r c 

(C7) 

or 

r 2 r 2 - r 2 t 2 

r r b _ r a r c 

(C8) 



ro h tr* 



t ro h tr* 
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Then, if equations (C4), (C 5 ), and (C 6 ) are substituted into equa- 
tion (C 8 ) there results 


lm(l - 2x)(m - 2y) = 0 (C 9 ) 

Infinite error therefore occurs if 


x = - (CIO) 

2 

or 

y = | (cn) 

Thus, infinite error occurs on the planes which are perpendicular 
bisectors of the sides of the rectangle. 
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APPENDIX D 

STATION-LOCATION ERRORS 


In order to determine the errors in position due to errors in sta- 
tion location, the following range differences, which follow the defini- 
tions, may be examined: 

s a = \/(x - x a ) 2 + (y - y a ) 2 + (z - z a ) 2 - '/x 2 + y 2 + z 2 (Dl) 

s b = \/(x - x b ) 2 + (y - y b ) 2 + (z - z b ) 2 - \/x 2 + y 2 + z 2 (D2) 

s c = \/(x - x c ) 2 + (y - y c ) 2 + (z - z c ) 2 - ^|x 2 + y 2 + z 2 (DJ) 

These equations may be interpreted as follows: s a , s^, and s c 

are measured input quantities; x a , y a , and so forth are constants 
depending on station locations; x, y, and z are the unknowns to be 
found. 

If it is assumed that exact data s a , s-^, and s c are being 
received, x, y, and z could be computed exactly if x a , y a , z a ; 
x b , y^, z b ; * • • are kno™ exactly. But if there are errors in 
measuring x a , y a , and so forth the problem is then to find the corre- 
sponding errors in x, y, and z. The procedure is to find the varia- 
tions in x, y, and z due to variations in x a , y a , and so forth 

with s a , and s c held constant. Taking the total differential 

of s a in equation (Dl) and setting it equal to zero gives 

ds a = 0 


ds s 

6xc 


dx a + 




6s„ 6s a 

d X a + T 1 dz a + + 

dz a ox 


Ss 

dy 


a 


dy 


OSr- 

c 

dz 


dz 




But, 


^ s a 

dx a 


- x„ 


(D5) 


X 


PO H 1 7 * 


ro h f 
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ds a = y - y a 


ds a _ z ~ z t 


Sz c 


Substituting equations (D5), (d6), and (D7) and a x , By, and a z 

ds a c>So ds a , . . 

- — , and into equation (D4) gives 

ox oy oz 


X - Xc 


y - y a , z - z e 


a x dx + ay dy + a z dz = — dx a + — dy a + 

** d Q *■ O 


dz c 


r a “ a ' ^a a 


Similarly, 


x - xv, y - y b z - z b 

b x dx + by dy + b z dz = — — - dx b + — — — dy b + — — — dz b 


x - x c . y - y c z - z c 


c Y dx + c v dy + c~ dz = — dx,. + ^ — dy_ + 2. dz r 

x y J 2 rc c r c r c 


These equations may be interpreted as follows: dx, dy, and dz 
the errors in x, y, and z, respectively, caused by the errors 
dy a , and so forth. 

Solving equations (D8), (D9), and (DIO) for dx yields 



+ 


z - z. 


**c 


dz c C x 


(d6) 

(D7) 

for 

(D8) 

(D9) 

(DIO) 

are 

te a , 


(Dll) 


28 


The probable errors in x a , and so forth may be expected to be 

equal but uncorrelated. Thus, letting aa be the probable error in 
any one of the quantities yields 


crx a = ax b = crx c = ay a = oy b = oy c = a z a = crz b = a z c = a a (D12) 


The probable error in x is the square root of the sum of the squares 
of each error term. Thus, 


• MV • [(Hr*) ! • (^r)’ * (^)>' • p*)’ * r)' * 


C 2 (D15) 


Substituting equations (10), (ll), and (12) into equation (D15) gives 


ox - f 0>lM 

In a similar manner, 

oy = 2S ^Ay 2 + By 2 + Cy 2 (D15) 

02 = ^ ^A Z 2 + B z 2 + C z 2 (Dl6) 

Since 

err = V^x^ 4- cfy^ + erz^ (D17) 

equations (Dl4), (D15)> and (Dl6) combine to become 


g = £ ^A x 2 + Bx 2 + C x 2 + Ay 2 + By 2 + C y 2 + A z 2 + B z 2 + C z 2 (Dl8) 

Since the right side of equation (Dl8) is identical to the right side of 
equation ( 38 ), 


H <M 



ro h t" 1 
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APPENDIX E 

DETERMINATION OF VELOCITY 


In order to find expressions for the velocity, consider deriva- 
tives with respect to time instead of differentials in equations ( 25 ), 
( 26 ), and ( 27 ) as follows: 


ds a 

dt 


(El) 


ds b 

dt~ 



dt 




dz 

dt 


(E2) 


ds c 


dt 





(E5) 


Solving for the x component of velocity gives 


dx = 1 
dt E 





(El) 


However, differentiation of equation ( 7 ) gives 


ds a 

dt 


- ( f ° 


fa)A 


(E5) 


Thus, 


Vx = dt = Et Ax(f ° “ fa) + Bx(f ° ' fb) + Cx(f ° " ^ 


(E6) 


The other components may be found in a similar manner. 
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APPENDIX F 
ERRORS IN VELOCITY 


The error in v x caused by errors in frequency measurement is 
obtained from equation (E6) (if it is assumed that position is known 
accurately) as follows: 


dv x = |[A X d(f Q - f a ) + B x d(f Q - f b ) + C x d(f D - f c )] (FI) 

Nov if it is assumed that the probable errors a(f Q - f a ) , 
d(f 0 “ f b)> and 0 ( f o “ fc) are all independent and equal to af , 
the equation may be written 


v x 



+ B x + c, 


(F2) 


Similarly, 




y e 


2 + By 2 + C y 2 


v z = | af\|A z 2 + B z 2 + C z ‘ 


The probable error in v* is then 


(F3) 

(F4) 


av = Mf M x ^ + B x " + Cx* - + A^ 


+ B + C 

y y 


+ a z + b z 


+ c. 


(F5) 


Thus, 


av = Aeaf 


^ H W r 



a H CVJ KN-4- 
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A V A 


i Infinite error occurs 
in these vertical planes 









(a) Altitude, 50*000 feet. 



(b) Altitude, 150,000 feet. 




(c) Altitude, 500,000 feet. (d) Altitude, 1,000,000 feet. 


Figure 8.- Error-factor contour plots for isosceles-triangle configura^ 
tion. Radii of charts, 1 x 10^ feet. 
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(a) Altitude, 50,000 feet. 


(b) Altitude, 150,000 feet, 




(c) Altitude, 500,000 feet. 


(d) Altitude, 1,000,000 feet. 


Figure 9.- Error-factor contour plots for particular configuration. 
Radii of charts, 1 X 10^ feet. Coordinates of stations in 
1,000 feet (origin at Wallops Island, Va.): station A, Dover, Del 
(-80, 454, -5); station B, Langley Field, Va. (-207, -310, -3); 
station C, Hatteras, N.C. (97 > -940, -21). 



(a) Altitude, 50,000 feet. 


(b) Altitude, 150,000 feet. 



(c) Altitude, 500,000 feet. (d) Altitude, 1,000,000 feet. 

Figure 10.- Propagation-error-factor contour plots for equilateral- 
triangle configuration. Radii of charts, 1 x 10^ feet. 


NASA - Langley Field, Va. L-12J4- 




